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Abstract 

We compute the quark-antiquark potential in three dimensional mas¬ 
sive Quantum Electrodynamics for arbitrary fermion mass. The result 
indicates that screening prevails for any quark masses, contrary to the 
classical expectations, generalizing our previous result obtained for large 
masses. We also test the validity of several approximation schemes using 
a detailed numerical analysis. The classical result is still reproduced for 
small separation of the quarks. 
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1 Introduction 


A proper study of the problem of screening and confinement is of consid¬ 
erable importance in our understanding of gauge theories. To avoid the 
complexities of four dimensions these studies are usually confined to lower 
dimensions. In this framework, a deep physical interpretation has been 
achieved. Indeed, in two-dimensional QEDQ, one obtains screening in 
the massless case, but confinement in the massive quark case, realizing the 
expected picture. 

For QCD in two dimensions Gross et al were the first to discuss 
the subject. If dynamical fermions and test charges are in different rep¬ 
resentations, they find screening or confinement in some particular cases 
depending on whether the fermion is massless or massive. A similar conclu¬ 
sion in an identical setting has been arrived at for the massless case in |^. 
If, on the other hand, all fermions are in the fundamental representation, 
then screening prevails independently of the quark mass @. 

General inquires in two dimensional gauge theories have been performed 
recently by several authors |^], concerning the 6 vacuum structure, screen¬ 
ing, confinement and chiral condensates. In three dimensions related ques¬ 
tions were studied in 

In three dimensional space-time, for an abelian gauge group, the ques¬ 
tion of screening versus confinement has been recently analysed for large 
fermion masses |0] in which case the fermionic determinant can be com¬ 
puted as a series in the inverse mass. The conclusion was that, contrary to 
classical expectations, the theory is in the screening phase. Although this 
is expected from the fact that a Chern-Simons term develops and there is 
a topological mass generation, it is a further indication that the dynamics 
of gauge fields and the deep problem of screening versus confinement is far 
from being settled by a simple inspection of the classical behaviour of the 
theory. In the case of three dimensional QED, the outcome reveals that 
the vacuum polarization is once more capable of developing configurations 
that screen the external quarks, presumably modifying the dynamics of 
quark-antiquark bound-states. 

Here we extend the analysis of our previous work |Q] in order to in¬ 
clude all values of the fermion mass parameter. An explicit expression 
for the quark-antiquark potential is obtained following the usual ideas of 
bosonisation|]^, |8| but an analytic form cannot be obtained, and we re¬ 
sort to the use of numerical analysis. The results show that the screening 
phase obtained in the large mass limit persists for any value of the mass 
parameter (including vanishing mass). Next, the validity of certain ap¬ 
proximation schemesQ is tested. Using these approximations a simple 
form of the quark-antiquark potential can be given, which is compared 
with the exact form. We still obtain screening; moreover the behaviour of 
the functions is very accurately described by the proposed approximations. 

The paper is divided as follows. In section two we present the computa- 
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tion of the potential and the numerical results. We draw the potential for 
different values of the mass parameter, showing that its form is essentially 
the same in the whole range of real values for the mass term, leading to the 
screening phase. Furthermore, we test the approximation forms of all the 
functions necessary for computing the potential. Section three is reserved 
for conclusions and discussions. 


2 Computation of the quark-anti quark po¬ 
tential 


The partition function of three dimensional massive QED in the covariant 
gauge, in the presence of an external source J^, is given by 


Z = I d[(l:,'il:^A^]6{d^A^)exp< i / — m — eyl)V’ 




( 1 ) 


where is the field tensor, = d^A^ — dyA^j,. 

The bosonised version of the above defined action in the weak coupling 
approximation is given by the expression []l^, [Tl|| 


Z 


J dA^5{dfj,A^)exp J d^xx 




( 2 ) 


where the dots stand for non quadratic terms in the gauge field Afj,. This 
result will describe the partition function of the Maxwell-Chern-Simons [I^] 
theory in the covariant gauge in the infinite mass limit 0, as we see from 
the explicit expression for the self-energy of the gauge field, 11^,^, given by 
the expression 


pH 

= H{p)iefj,yp^ + 


G{p) +p^ 


9ijlv 


p2 


where the functions G and H are given by the forms 


(3) 


H{p) 

G{p) 


e^p^ 

dvr 

-P^ - 


dt- 


e^p^ 


[m^ 


dt 


t{l — t)p‘^YG ’ 


[m^ — — f)p2]i/2 


(4) 

(5) 


and p = \/—p^- We compute the potential as being the difference be¬ 
tween the Hamiltonian with and without a pair of static external charges 
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separated by a distance L, 


V{L) = Hg - Ho =-{Lg - Lo) 

= —q J cfxAfj_5^^[6{x^ + L/2)S{x‘^) — 6{x^ — L/2)6{x'^)] 

= -q[Ao{x^ = -L/2, = 0) - Ao{x^ = L/2, = 0)] , (6) 


where we have integrated over the two space components in order to find 
the potential, and considered the source as corresponding to two fixed 
charges of magnitude q located at the points defined by the respective 
delta functions. Note that Lq(Lo) denote the Lagrangeans in the presence 
(absence) of the charges. 

We now consider the equations of motion associated with the action 
defined by means of (|2|). The field equation in the covariant gauge reads 

- H{d = + G(9) + J;. = 0 (7) 

Defining the curl of as 

( 8 ) 


the equation of motion can be expressed as 


[^nonloc T ^nonloc\Afj, — J fnonlocJfi 


where 


^nonloc 

^nonloc 


fnonloc 


G{d) 

H^{d) 

d‘^G{d) 

H{d) 

G{d) 


(9) 


( 10 ) 


In the absence of sources, and in the large m limit Q, it reproduces the 
familiar massive mode of Maxwell-Chern-Simons theory [p)2| . From (^) it is 
seen that an expression for Aq is required to calculate the potential. This 
is given in terms of the curl (|^ by 


A 2 = —diAo 


( 11 ) 


The time independent solution for A 2 corresponding to the sources 
describing static quarks can be obtained from (|^. Using this result with 
( [III) finally yields, after integrating over the angular variables, 

Mt,L) = A„{a,L) = ArG^M2^ik ( 12 ) 
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The integration over the angular variables in the Fourier transformation 
led to the Bessel function Jo^kL); in the case where the denominator is 
given by the familiar result (i.e. Feynman propagator, which also appears 
in the large mass limit, see |^) the result of the integration is just the 
modified Bessel function 

The potential is now found from (^), (|^ and (p^), reading 


V{L) 


q 2 .00 

vr Jo 


kJo{kL) 


dk 


(13) 


where the functions G{k) and H{k) are given by the expressions 


G{k) 

H{k) 


r + 


e^m 


dvr 


/ 2m 

k \ 

J k \ 

1 _ - 

— 

arctan - 

L V k 

2m / 

V2myJ 


e^mk 

‘It: 


arctan 


k 

2m 


(14) 

(15) 


The above equation takes a particularly simple form in the infinite mass 
limit, 


ViL) = — 



^1 + ^ 


(16) 


This reproduces our earlier results in Q. The asymptotic form of the 
Bessel function signalizes screening. 

For arbitrary mass however, a simple closed form expression cannot be 
obtained. We therefore have to use numerical methods. They are presented 
as follows. We first plot the function V{L) given in (1^) as a function of 
L for various values of the mass parameter m. The result is plotted in 
figure 1. It is immediately obvious that the screening effect is qualitatively 
independent of the mass, and the quantitative dependence extremely small. 
Indeed, the graphs are almost bound inside a rather narrow band defined 
by the results obtained for m = 0 and m = 00 . 

In general, due to the appearance of non algebraic functions the ex¬ 
pressions appearing in (14) and (15) are rather clumsy. In reference |^] 
simple expressions have been derived which, according to the authors, give 
a good approximation to these functions in the whole range of values of 
the parameter. The approximations are 


G{k) 

H{k) 






1 

2 



e^mk"^ 




1 

2 


(17) 

(18) 


We tested this assumption for the computation of the potential, com¬ 
paring the asymptotic result with the one obtained with the approxima¬ 
tions for m = 1000. We repeated the procedure for m = 0.1, which shows 
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Figure 1: Potential as a function of the distance for various values of the mass param¬ 
eter. The value m = 1 practieally coincides with the asymptotic value m = oo. 


similar findings. The result is shown in figure |2|, indicating that the approx¬ 
imation agrees remarkably well with the exact results. We also checked the 
approximations directly in figure Q. Comparison of the expressions for G e 
H using the approximations and the exact result shows that there is little 
discrepancy. 

The approach to asymptotics in the computation of the potential has 
also been analysed. We have verified that it is very quick. Indeed, for rea¬ 
sonably low values of the mass the potential already shows the asymptotic 
value. We illustrated the behaviour in the case m = 10 in figure which 
was done with the exact expressions. 


3 Discussion 

Here we worked out an approach for obtainning the semiclassical inter¬ 
quark potential for arbitrary values of the fermion mass parameter, gen¬ 
eralizing our previous workQ, where only the infinite mass approximation 
was analysed. In that case, the expression for the effective action is local, 
and the interquark potential could be computed in closed form, showing 
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Figure 2: (Left) Potential function for infinite mass, and for m = 1000 using and 

(Right) Potential function for m = 0.1 using the exact results & and and 

approximate forms and 



Figure 3: Comparison between the potential for m = 10 and the asymptotic result. 


explicitly that the model lives in a screening phase. We attempted here to 
go beyond the large mass limit. However, the expressions thus obtained 
are nonlocal, and we had to resort to numerical simulations. The results 
are nevertheless rewarding, especially in view of the extremely mild depen¬ 
dency upon the fermion mass, i.e. almost every physical quantity related to 
the screening potential is almost independent of the mass, for 0 < m < 1 
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Figure 4: Function G (left) and H (right) with m = I, using the exact definitions & 
and and the approximations 0) and &■ 


(we suppose e = 1) and achieves the asymptotic value already for m of 
order unit. 

The screening obtained for all values of the mass parameter supports 
the observations obtained in two-dimensional QCD|Q, where the screening 
phase also prevails almost universally (see also ||^). This opens up the 
discussion for a large number of interesting possibilities. In particular, it 
is interesting to stress that the same mechanism may work for the non- 
abelian case in three dimensions, since in the large mass limit the effective 
action turns out to be the non abelian generalization of the Maxwell Chern 
Simons theory. In an axial gauge, and in the weak coupling limit, the 
Maxwell Chern Simons action coincides with the abelian counterpart, and 
the same conclusions are expected. If one dares speculate that the mass 
dependence is as mild as we have obtained in the above discussion, then all 
conclusions can be carried over to the non abelian case as well, a tantalising 
result! This would imply an almost universal screening behaviour for low 
dimensional systems. We hope to come back to these interesting questions 
in a future work. 

We also tested the approximative formulae (|T7|) and (jT^) , usually taken 
as a good approximation within 10% accuracy. The difference between the 
potential calculated using this approximations and the exact results is less 
than 10“^, too small to be seen in figure 

We finally comment on the infinite and zero mass limits of the model. 
In the latter case we have to go beyond the quadratic approximation, 
since higher corrections have to be computed in a gauge theory where 
further powers of the external momenta show up in the computation of 
diagrams. One loop fermionic diagrams in gauge theories result in powers 












of momentum and functions of momentum over mass, in general /(^). 
While the limit j? ^ 0 and m —>■ oo is unambiguous, the double limit 
p —> 0 and m ^ 0 is not well defined, and depends on the order they are 
taken. Therefore strong infrared dependence on the mass may invalidate 
the procedure. Notice that the discussion of screening pressuposes a large 
distance (namely small momentum) limit, which may not commute with 
the zero mass limit. The infinite mass limit obtained in the quadratic 
approximation is however expected to survive even in the non-quadratic 
regime. We hope to come back to these points in a future publication. 
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